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Abstract 

A concise and critical review of B aimer-dominated shocks (BDSs) is presented, summarizing the state of theory and observations, 
including models with/without shock precursors and their synergy with atomic physics. Observations of BDSs in supernova rem- 
nants are reviewed on an object-by-object basis. The relevance of BDSs towards understanding the acceleration of cosmic rays in 
shocks is emphasized. Probable and possible detections of BDSs in astrophysical objects other than supernova remnants, including 
pulsar wind nebulae and high-redshift galaxies, are described. The case for the continued future of studying BDSs in astrophysics 
is made, including their relevance towards understanding electron-ion temperature equilibration in collisionless shocks. 

Key words: acceleration of particles - atomic processes - cosmic rays - galaxies: high-redshift - hydrodynamics - ISM: general 
- line: formation - radiation mechanisms: general - shock waves - stars: winds, outflows 



1. Introduction 



Shocks are present wherever matter is accelerated past the 
sound speed of the medium they are propagating through 
Their ubiquity precludes any one author from 
performing an exhaustive review of their study, which includes 
theory, observation, experiment and simulation — both terres- 
trial and astrophysical. In this review, we focus on a specific 
class of astrophysical shocks known as "B aimer-dominated 
shocks" (BDSs; Figure [U, previously reviewed by Raymond 
SB Traditionally observed around historical supernova 
remnants (SNRs), these are ~ 200-9000 km s" 1 shocks imping- 
ing upon the ~ 0.1-1 cm -3 interstellar medium (ISM), mani- 
festing themselves as faint (but limb-brightened), optical fila- 
ments. 

BDSs are characterized by: 

• Having relatively high velocities (> 200 km s" 1 ) such that 
they are "non-radiative" (©; 

• Presence of strong hydrogen lines with narrow (~ 10 km 
s" 1 ) and broad (~ 1000 km s" 1 ) components (e.g., Figure 

0; 

• Absence or weak presence of forbidden lines of lowly- 
ionized metals; 

• General lack of non-thermal X-ray emission at the loca- 
tions where two-component Ha lines are detected (but see 

That the forbidden metal lines are weak relative to the hydrogen 
(and helium) lines is simply a consequence of the high temper- 
atures in BDSs allowing the production of both > 10 eV and 
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Figure 1: Schematic structure of a B aimer-dominated shock. As the cold, pre- 
shock hydrogen atoms stream past the collisionless shock front, they interact 
with the post-shock protons via charge transfer reactions to create a secondary 
population of atoms known as "broad neutrals" ( 36.31 . The post-shock elec- 
trons and protons are heated in a thin layer with thickness ~ / gyro (equation 
fill ). Broad neutrals carry information about the post-shock gas and are cre- 
ated in the post-shock region known as the "shock transition zone" ( 36.lt with 
a thickness ~ /atomic (equation flOl ). The densities of the pre- shock hydrogen 
atoms, broad neutrals and protons are shown; the density jump for the protons 
reflects the Rankine-Hugoniot conditions. The dashed arrows indicate that cos- 
mic rays and broad neutrals may diffuse upstream and heat the pre-shock gas 
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Figure 2: Example of the two-component Ha line from a B aimer-dominated 
shock located in northeastern RCW 86, a Galactic supernova remnant. The ver- 
tical arrows indicate [N n] i6548, 6584 emission. Courtesy of Eveline Helder. 



< 1 eV lines — the relative line strengths then scale as the el- 
emental abundance — in contrast to the conditions in emission 
line nebulae where the hydrogen lines are suppressed because 
of low temperatures. BDSs are somewhat simple to interpret, 
because radiative cooling and recombination are unimportant in 
the immediate post-shock region. 

Combined with proper motion measurements, BDSs can be 
used to infer distance^ to spatially-resolved objects. If both 
the narrow and broad components are detected, the ratio of 
electron-to-ion temp eratures just behind the shock front can 
be inferred l34ill0lll . thereby providing constraints for plasma 
physics simulations of electron-ion^ equilibration. BDSs can 
be used to extract information on cosmic ray ions and study 
particle acceleration in partially neutral media. The broad hy- 
drogen lines produced generally serve as direct tracers of the 
post-shock proton distribution. Understanding the physics of 
these fast, spatially-resolved shocks allows us to better under- 
stand if shocks are important in distant, spatially-unresolved, 
astrophysical objects. 

In §|2l we develop some intuition for the physical conditions 
in BDSs by familiarizing the reader with order-of-magnitude 
time and length scales. In ^3] we gain historical perspective by 
tracking the pioneering efforts to study BDSs in SNRs. In 
we compile and summarize the SNR observations on an object- 
by-object basis. The relevance of BDSs towards studying cos- 
mic ray acceleration is emphasized in §5\ In we review the 
theoretical models of BDSs, including those with and without 
shock precursors; we also describe the atomic cross sections 
that are imperative for constructing these models. Probable and 
possible observations of BDSs in other astrophysical objects are 
described in ^7) The case for the continued future of studying 
BDSs is made in ^8) Table Q] contains a list of symbols that are 
commonly used throughout the paper. 



*It can be argued that, at least for SNRs, this method is an improvement 
over using the H-D relation to obtain distances 

2 In this paper, we use "ions" to mean non-leptonic ions, i.e., excluding elec- 



2. Order-of-Magnitude Physics: Time and Length Scales 

A shock with a velocity impacting upon ambient ISM with 
a density n, magnetic field strength T^sm and average gas mass 
yu has sonic and AlfvenicH Mach numbers of 



M s 



~ 100 



1000 km s 



- l )\m n ) U0 4 K/ 



and 



500 



1000 km s~ 



1 cm -3 m\\ 



1/2 



-1/2 



(1) 



(2) 



implying that high- velocity shocks are inaccessible to terrestrial 
investigations, where mu is the mass of the hydrogen atom and 
k is the Boltzmann constant. 

If the pre- shock hydrogen gas has a temperature To and emits 
a thermal, Gaussian-shaped line, its full width at half-maximum 
(FWHM) is 



W = 2 



2\n2kT { 



o 



21 kms 



m H 



\10 4 K/ 



(3) 



Since To ~ 10 4 is the temperature above which collisional ion- 
ization overwhelms recombination (e.g., see Figure 12.4 of Os- 
terbrock & Ferland 1 67]), narrow hydrogen lines with widths 
greater than about 20 km s" 1 are believed to be non-thermally 
broadened (see 35.11) . In this paper, we take the thermal veloc- 
ity of an ion to be the root mean square (rms) velocity in three 
dimensions, given by 



Vth,i 



(4) 



where T t and m; are the ion temperature and mass, respectively. 
The thermal velocity of an electron is obtained by inserting T e 
and m e into equation ©. 

The shock heats up the impacted material to a temperature of 



(y+\f\ki s \m p 



1000 kms" 



(5) 

where y is the adiabatic index of the shocked gas and m p is 
the proton mass. The second equality assumes y = 5/3, i.e., a 
monoatomic, non-relativistic gas. Emission lines from the post- 
shock gas typically have a width of 



W 



4v 



^Vin2 (r-D. 



(6) 



For hydrogen lines (see ®, there is a non-negligible depen- 
dence of W on 13 = T e /T p , the ratio of electron-to-proton tem- 
peratures EHfioi]. For y = 5/3, we have W ~ v s . 



3 We differentiate between Alfvenic and magnetosonic Mach numbers; the 
latter is defined as v s relative to the sound and Alfven speeds summed in quadra- 
ture. 
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Symbol 


Meaning 




shock velocity 


Av 


relative velocity between atoms and ions 


n 


ambient density (total) 


/ion 


pre- shock ion fraction 


^ISM 


magnetic field strength of ambient ISM 


To 

1 u 


pre- shock temperature 


- 1 s 


post- shock temperature 


T e , T p , Ti 


post-shock electron, proton and ion temperatures 


p 


ratio of electron to proton temperatures 


Wo 


full width at half-maximum of narrow hydrogen^ line 


W 


full width at half-maximum of broad hydrogen^ line 


hUn 


ratio of broad to narrow hydrogen^ line intensities 



f : typically Ha unless otherwise noted. 



Table 1 : Commonly-used symbols 



For lack of a better term, BDSs are described as being "non- 
radiative", meaning that the age of the shock is less than both 
the cooling and recombination time scales in the intermediate 
post-shock region, 



'age 



< minUcooi,^ 



(7) 



In other words, the shock does not produce radiation that is 
dynamically important. 

Accurate determinations of the cooling time scale requires 
specifying the elemental abundance of — and calculating the 
radiative transfer of photons through — the post-shock plasma 
, but a simple estimate can be obtained using 



^cool 



nA 



4xl0 5 yrs 

y \10 7 K/\lcm 



cm -3 10 -22 erg cm 3 s~ 



-l 



(8) 

where T s is the post-shock temperature, n is the ambient density 
and A is the cooling function. 

Recombination becomes important when the post-shock 
plasma cools to ~ 10 4 K; the recombination time scale is ap- 
proximately 



1 



na n 



8 x io 4 yrs ( ( -^-f 1 , (9) 
y VI cm" 3 / \10 4 K/ W 



where a rec is the recombination coeffic ient II47L 16711 . More ac- 
curate expressions for a YQC are given in 

It is apparent that for shocks with ages ~ 100-1000 years 
and velocities ~ 1000 km s" 1 impinging upon the ISM with 
~ 0.1-1 cm -3 densities, the non-radiative condition described 
in equation © is trivially fulfilled. 

BDSs are also collisionless shocks, because the length scale 
for atomic interactions at ~ 1000 km s" 1 (typically dominated 
by charge transfer reactions at these velocities; see 1 



/el 



1 



^^"atomic 



io 15 



cm 



(lcm- 3 ) 



(10) 



is much larger than reasonable estimates for the proton gyrora- 
dius, 



/gyro ~ 10 10 cm 



1000 km s" 



#ISM 
lyUG 



(ID 



where v p is the velocity of the proton. The corresponding gy- 
roradius for an electron is smaller by a factor of m e /m p , where 
m e is the electron mass. 

The condition / gyro <*c /atomic implies that the relative tem- 
peratures of the electrons and ions are determined by collective 
electromagnetic interactions, rather than by physical collisions. 
The time scale for Coulomb equilibration II971 l27ll is 



c 



cm" 3 10 



(12) 



where T e is the electron temperature and C is the Coulomb 
logarithm. Since t e { » t agQ , Coulomb interactions alone can- 
not bring about temperature equilibration in young SNRs (with 
ages ~ 10 3 yrs). The range of electron-to-ion temperature ratios 
is expected to be m e /m[ < T e /Ti < 1 II 1 1 Oil . 

In the absence of limb brightening, the intensity or brightness 
of a shock filament can be used to provide an estimate for the 
ambient density, since 



ha ~ 



HV S €ug Eii a 

4n 



= 5 x 10 6 erg cm 1 s 1 sr 



-2 -1 -1 



1 cm" 3 1000 km s" 1 02)' 
(13) 



where €n a is the number of Ha photons per ionization event and 
Etta ~ 2 eV is the energy of a Ha photon. 

Energetic particles and/or photons may be generated in the 
immediate post-shock region of a non-radiative shock. When 
the typical propagation velocities of the energetic particles ex- 
ceed the shock velocity — a condition automatically satisfied 
by the photons — they will travel upstream into the pre-shock 
gas and heat it, thus broadening the line width Wo. Such broad- 
ening mechanisms are known as "precursors" and may be re- 
sponsible for the unusually large values of Wo observed in sev- 
eral SNRs (see 35.11) . In shocks where a cosmic ray precursor 



4 The pre-shock hydrogen atoms affected by the precursor must still cross 
the collisionless shock front to produce the observed narrow lines. 
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is believed to be present, the length scale for cosmic ray prop- 
agation (/cr ~ kcr/v s ) must be shorter than that for ionization, 
thereby setting an upper limit to the cosmic ray diffusion coef- 
ficient: 

*cr < 10 24 cm 2 s" 1 (r^Y (— £ r) , (14) 

VI cm" 3 / \ 1000 km s" 1 / 

where the typical cross section for ionization is taken to be ~ 
10" 16 cm 2 . 

3. Historical Developments: the Remnants of Tycho, SN 
1006 and the Cygnus Loop 

In 1959, Minkowski Il66ll announced, at the Paris symposium 
of the International Astronomical Union, the detection of bright 
Ha filaments around the remnants of Kepler, Tycho and the 
Cygnus LoopH. Following an unsuccessful search for the opti- 
cal remnant of supernova (SN) 1006 by Walter Baade, van den 
Bergh lll02h in 1976 reported what he termed "delicate wisps 
of filamentary nebulosity" about 10' from the center of the as- 
sociated radio source. He noticed that the morphology of the 
filaments were similar to those detected in Tycho 's SNR (the 
remnant of SN 1572) and S147. Shortly afterwards, Kirshner 
& Taylor 115011 presented Ha velocity measurements of up to 
300 km s" 1 in the Cygnus Loop equivalent to shock tempera- 
tures ~ 10 6 K, consistent with the observation of soft X-rays 
and [Fe xiv] /15303 emission behind the shock front. 

In 1978, McKee, Cowie & Ostriker Q suggested that the 
high-velocity Ha emission may be from pre-existing clouds in 
the SNR being accelerated and shock-heated by the blast waves, 
although they leave open the question of whether the clouds can 
survive in the environment of the SNR. One prediction of the 
"accelerated clouds" model is that the Ha emission should be 
accompanied by substantial emission in, for example, the [O i] 
,16300 and [O m] ,15007 lines. 

In the same year, an optical spectrum of Tycho 's SNR ob- 
tained by Kirshner & Chevalier Dill provided important clues 
towards understanding the emission mechanisms associated 
with the filaments. Modeling the thermal X-ray continuum 
measured by Ariel 5 yielded emission components with two 
different temperatures, 0.55 and 3.5 keV. Assuming a distance 
to Tycho 's SNR of 6 kpc, the tangential velocity of the fila- 
ments was estimated to be 5600 km s" 1 , corresponding to post- 
shock temperatures near 40 keV. A natural way to reconcile 
these two observables was for incomplete temperature equilib- 
rium between the electrons and ions to occur — in this case, for 
Ti/T e ~ 20-30 such that the electron temperature is T e ~ 5x 10 5 
K. 

A second important clue provided by the observations of Kir- 
shner & Chevalier was that only Ha and Hj3 emission lines 
were observed over the wavelength range of 4700-8000 A. 
Forbidden lines of singly- or doubly-ionized metals (e.g., [O 
m] /14959, 5007) were conspicuously absent, thereby preclud- 
ing their use as density or temperature diagnostics and placing 



5 Minkowski also noted the presence of more than 200 "small condensa- 
tions" in the nebulosity around Cassiopeia A. 
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Figure 3: Diagnostics for edge-on shoc ks using the Ha line, based on calcu- 
lations from van Adelsberg et al. llOlil . The width of the broad line W and 
the ratio of the broad to narrow line intensities 4 //„ both depend on the shock 
velocity v s and the ratio of electron to proton temperatures p. For illustration, 
we adopt Case A and B conditions for the broad and narrow lines (see 36.3b , 
respectively, and fix the pre-shock ion fraction at f[ on = 0.5. 



into question the accelerated clouds model. Two interpretations 
were suggested: either the forbidden lines were collisionally 
quenched due to high densities, or the high temperatures caused 
all of the metals to be in their highest ionization states. The ab- 
sence of the [O m] ,15007 line and the [O n] ,17319, 7329 dou- 
blet requires densities in excess of 10 7 cm -3 , orders of mag- 
nitude above typical ISM densities, thus ruling out the former 
possibility. 

Also in 1978, Schweizer & Lasker l87ll obtained spectra in 
the wavelength range of 3650-7500 A for BDSs at the north- 
western rim of the remnant of SN 1006. Only B aimer emission 
lines were detected; the [O i], [O n], [O m], [N n] and [S n] lines 
were either absent or weak (< 8% by intensity relative to Ha 
when uncorrected for reddening). 

An alternative model was proposed by Chevalier & Raymond 
W18I1 . in which a collisionless shock impacts partially neutral 
gas. Their model was able to explain the intensity, spectrum 
and width of the filaments observed in Tycho 's SNR, including 
the weakness of the forbidden metal lines. Lines of neutral he- 
lium were predicted to be the next brightest optical lines after 
the hydrogen Balmer lines; for example, the He i ,15876 line is 
expected to be about 10 times fainter than the HJ3 line. Limb 
brightening increased the intensity at the edge of the SNR by a 
factor ~ C^snr Aon) 1/2 , where 7?snr is the radius of the remnant 
and /ion ~ 10 16 (1 cm -3 /ft) cm is the ionization length scale. By- 
chkov & Lebedev fl3ll estimated the effect of photoionization 
of the pre-shock gas by photons from the post- shock gas; they 
predicted that high-velocity hydrogen emission should be ex- 
pected from Kepler's SNR, but not from nebulae surrounded by 
H ii regions (e.g., Cas A, Crab). 

In a seminal 1980 paper, Chevalier, Kirshner & Raymond 
1 19] reported the detection of Ha, Hfi and Hy lines in Tycho's 
SNR and also expanded upon the shock wave theory proposed 
by Chevalier & Raymond. The Ha line appeared to have two 
components: an unresolved, narrow component due to direct 
excitation of cold, pre-shock hydrogen atoms; and a broad com- 
ponent with a FWHM of W = 1800 ± 200 kms" 1 . The broad 
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line component arises from a secondary population of hydro- 
gen atoms created from charge transfer reactions between cold, 
pre-shock atoms and hot, post-shock protons that are called 
"broad neutrals". Thus, broad neutrals carry valuable infor- 
mation about the post-shock gas. 

The Chevalier, Kirshner & Raymond paper established sev- 
eral important results. Firstly, W is a diagnostic for the shock 
velocity (equation (6)). Secondly, the ratio of broad-to-narrow 
line intensities, denoted by 4/4, is a sensitive function of the 
shock velocity, because the charge transfer cross section de- 
pends strongly on the electron and proton velocities (see 36.31) . 
Both W and 4/4 depend on the degree of temperature equi- 
libration between electrons and ions (see Figure 0. Thirdly, 
it was realized that the conversion of Lyman to B aimer photons 
via scattering, known as "Lyman line trapping" (see 36.3l for de- 
tails), is a non-negligible effect in Tycho's SNR. In pre-shock 
media where the optical depth to Lyfi photons equals or exceeds 
unity, one Ha photon is produced for every 5-10 atoms that are 
ionized, a result that is confirmed by later calculations ll65. 411. 
Combining the inferred value of the shock velocity with proper 
motion measurements, the distance to Tycho's SNR was revised 
downward to 2.3 ± 0.5 kpc. 



4. Observations of Balmer-Dominated Shocks in Historical 
Supernova Remnants 

Following the pioneering studies described in §|3] the field 
of B aimer-dominated shocks grew out of optical observations 
made of faint shock filaments around historical SNR^l that pre- 
dominantly have Type la (i.e., thermonuclear) SN progenitors. 
A sample image of the shock filaments is shown in Figure |4] for 
RCW 86. Reported values of W and 4 /4 are typically obtained 
from fitting Gaussians to the observed lines. 

The translation of the measured W and 4 /4 values (usually 
from Ha lines) into reported values of v s is traditionally ham- 
pered by uncertainties in the temperature equilibration between 
electrons and ions. Generally, the geometry of the ambient 
magnetic field in these SNRs is poorly known, particularly in 
remnants where radio emission is weak. 

4.1. Galactic SNRs 

4.1.1. The Cy gnus Loop 

Raymond et al. J74ll report optical and ultraviolet observa- 
tions of shock filaments at the eastern limb of the Cygnus Loop. 
The two-component Ha line is detected; by fitting Gaussians 
to each component, the narrow and broad FWHMs are mea- 
sured to be 31 and 167 km s" 1 , respectively. The centers of the 
lines differ by less than the uncertainty in centering the broad 
component, indicating that the shock is moving almost exactly 
transverse to the line of sight, i.e., an edge-on shock. A search 
for the variation of the broad Ha line width across the shock 
filament proves to be unsuccessful. 



Raymond et al. construct two models: the first assumes that 
plasma turbulence equilibrates the electron and ion tempera- 
tures (i.e., T e = Ti) rapidly (compared to any time scale for 
atomic interactions); the second assumes only (slow) Coulomb 
equilibration, i.e., T e <*c 4. They demonstrate (see their Table 
2) that the models with only Coulomb equilibration match a se- 
ries of line intensities, as well as the strength of the hydrogen 
two-photon continuum, better than the rapid equlibration mod- 
els. Their models also predict that bright O vi emission should 
be present in ~ 200 km s" 1 shocks. 

Fesen & Itoh [25] measure the line intensities of H, He 1, 
He 11 and several low-ionization metals along a different portion 
of the same filament studied by Raymond et al. By matching 
the theoretical and observed line intensities, they conclude that 
two- and three-temperature models fit the data better than the 
one-temperature model. The three-temperature models are mo- 
tivated by wishing to distinguish between shocked electrons and 
secondary electrons produced via ionization of neutral atoms. 
Nevertheless, the three-temperature model does not provide 
significant improvements over the two-temperature one. Cu- 
riously, the observed [N 11] ^6583 and [S 11] ^16717,6731 line 
intensities exceed the model predictions by at least an order of 
magnitude, which motivate two interpretations: either the ob- 
served shock has a radiative component, or dust grains present 
in the non-radiative shock are reflecting emission from a proxi- 
mate, radiative filament. 

Ghavamian et al. 113 111 measure W = 262 ± 32 km s" 1 
and 4/4 = 0-59 ± 0.3 in northeastern Cygnus. Using their 
models (see 36.11) . they find nearly complete (J3 « 0.7-1) 
electron-proton temperature equilibration based on the HJ3 line 
(4/4 = 0-99 ± 0.3), but are not able to match the measured 
4/4 value for Ha. 

Salvesen, Raymond & Edgar l85ll use proper motion mea- 
surements and X-ray spectra to constrain the ratio of cosmic ray 
to gas pressure behind BDSs in the northeastern rim. Adopting 
a distance of 576 ± 61 pc to the Cygnus Loop, they derive shock 
velocities from the proper motion data. Two-temperature fits to 
ROSAT spectra yield post- shock temperatures; the lower of the 
two temperatures is used to obtain an upper limit to the cosmic 
ray pressure. Full equilibration between electrons and ions is 
assumedQ For 18 non-radiative Ha filaments, they find that the 
ratio of cosmic ray to gas pressure is close to zero. 

4.7.2. The Remnant of SN 1006 

More than six years following an unsuccessful search by 
B.M. Lasker, Kirshner, Winkler & Chevalier [52] report a con- 
vincing detection of the two-component Ha line from a fila- 
ment along the northwestern rim of SN 1006. They measure 
W = 2600 ± 100 km s" 1 and 4/4 = 0.77 ± 0.08. Brack- 
eting the inferred shock velocities by models with no and full 
temperature equilibration between electrons and ions, and com- 
bining them with proper motion measurements, the distance to 
SN 1006 is estimated to be in the range 1.4-2.1 kpc. Smith et 



6 We refrain from using the term "B aimer-dominated supernova remnants", 
because a multitude of both non-radiative and radiative shocks may exist around 
a given remnant. 



7 This is a conservative assumption for the purpose of obtaining upper limits 
to the cosmic ray pressure, which is dominated by the ions since they constitute 
the bulk of the fluid mass. 
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al. W93I1 later reanalyze the spectrum from Kirshner, Winkler & 
Chevalier using improved line fitting techniques, thereby refin- 
ing the broad Ha line width to 2310 ± 210 km s" 1 and 4/4 to 
0.73 ± 0.06. 

Laming et al. [58] report estimates of the ratio of electron 
to proton temperatures behind a SNR shock front, by analyzing 
ultraviolet lines of He n ,11640, C iv ,11550, N v ,11240 and O vi 
,11036, and concluding that complete electron-ion equilibration 
is disfavored (see also $8.11) . 

Ghavamian et al. [32] report the first detection of the He i 
,16678 line in a BDS, in addition to measuring the broad Ha 
line width (W = 2290 ± 80 km s" 1 ) and a possible (~ 1.5cr) 
detection of the He n ,14686 line. The broad to narrow line ratio 
is measured for both Ha (I b /I n = 0.84+°^) and HJ3 (I b /I n = 
0-93+o Based on their models (see lO . they infer J3 < 
0.07. ' 

Raymond et al. lF78h use a deep Hubble Space Telescope im- 
age to resolve the Ha filament in the northwest. At 40 positions 
along the filament, they extract spatial profiles of the Ha in- 
tensity across the shock front. The exact details of the profile 
depend upon a combination of the pre- shock neutral density 
and neutral fraction — thus, matching the profiles to models 
lf78l l42ll allows these pair of values to be determined at each 
position. Taking the neutral fraction to be 0.05-0.2, Raymond 
et al. infer the pre-shock neutral density to be about 0.15-0.4 
cm -3 (see also §5 of |42|]). Also detected are parsec-scale rip- 
ples in the filament, which may reflect ~ 20% density fluctua- 
tions caused by interstellar turbulence. 

The distance to SN 1006 is a subject of controversy with 
quoted values ranging from 0.7 to 2.5 kpc; these estimates are 
anchored by the fact that spectra of the Schweizer-Middleditch 
star 118 811 . located at ~ 1.05-2.1 kpc, exhibit absorption by SN 
1006 ejecta 1 12]. Hamilton, Fesen & Blair [37] measure the ex- 
pansion velocity of ejecta at the reverse shock to be 7026 ±10 
km s" 1 , implying a shock radius of about 7 pc. Since the rem- 
nant is approximately 15' across, this corresponds to a mini- 
mum distance of about 1.6 kpc. The most direct determination^ 
is to combine measurements of the broad Ha line with filamen- 
tary proper motions lll08ll . Van Adelsberg et al. JlOlll use their 
estimates of the shock velocity to revise the inferred distance 
from ~ 2.2 kpc to ~ 1.6 kpc, a somewhat more comfortable 
figure in comparison to that of the Schweizer-Middleditch star. 

Overall, the remnant of SN 1006 has the distinction of be- 
ing the only instance of a spatially-resolved SNR where both 
efficient and inefficient sites for particle acceleration have been 
identified in BDSs (see $5.31) . but yet shows a lack of evidence 
for precursor heating ( $5.11) . 

4.1.3. Tycho's Remnant (SNR 1572) 

Kirshner, Winkler & Chevalier [52] detect the two- 
component Ha line from Tycho's SNR, measuring W = 1800 ± 
100 km s" 1 and 4 jl n - 1 .08 ± 0.16, thereby providing improve- 
ments to the measurements of Chevalier, Kirshner & Raymond 
and an estimate for the distance to the remnant of 2.0-2.8 kpc. 



8 The caveat is that the main systematic uncertainty (~ 10-30%) is contained 
within the atomic cross sections (see $6.3K 



Smith et al. [93] report significant variations in the Ha 
line profile along the filament containing "knot g", an over- 
dense region located at the northeastern rim of the remnant. 
The broad line width measured (W = 1900 ± 300 km s" 1 ) 
is similar to that reported by Kirshner, Winkler & Chevalier, 
but the broad-to-narrow line intensity ratio is somewhat lower 
(hi I n = 0.77 ± 0.09). The lower limit of the distance to Ty- 
cho's remnant is revised to 1.5 kpc, with the full range (1.5-3.1 
kpc) again being bracketed by uncertainties in the electron-ion 
temperature equilibration. 

Ghavamian et al. [31] report W = 1765 ± 110,2105 ± 130 
km s" 1 and 4/4 = 0-67 ± 0.1,0.75 ± 0.1 from two locations 
in knot g. They are able to measure B aimer decrements (i.e., 
the Ha/H/3 intensity ratio; see 36.3.11) for both the narrow and 
broad lines, finding Ha/H/3 « 4-6 (narrow) and « 3-4 (broad). 
The larger values of the B aimer decrements for the narrow lines 
indicate that Lyman line trapping, i.e., the conversion of Ly/3 to 
Ha and Lyy to H/3, is relatively more important for the narrow 
compared to the broad lines. 

Lee et al. [ 59] present evidence for a photoionization precur- 
sor, with a width ~ 10 16 cm, in knot g — He 11 /1304 photons in 
the post-shock region diffuse upstream and heat the pre-shock 
gas. The narrow Ha line has an unusually broad width (~ 45 
km s" 1 ), but the photoionization precursor is not believed to be 
the main cause of the broadening (see 35.11) . 

4.1.4. Kepler's Remnant (SNR 1604) 

Fesen et al. J26ll discover BDSs along the northern limb of 
Kepler's remnant. Averaging over BDSs observed at two slit lo- 
cations, they find W = 1750+200 km s" 1 and 4/4 = 1.1+0.25. 
The diffuse nature of the northern filaments is in contrast to the 
sharply-defined filaments observed in the remnants of Tycho 
and SN 1006. 

Blair, Long & Vancura ^ scrutinize the remnant in greater 
detail and confirm the result of Fesen et al. — diffuse BDS 
emission is present in the north indicative of the SN shock 
encountering lower density, uniform, pre- supernova material. 
They also report the discovery of the first BDSs with knotty 
structures. In general, the knots are categorized to be radia- 
tive, non-radiative and transitional; the transitional knots have 
weaker [S 11] emission compared to the bright, radiative knots, 
while the best examples of the non-radiative knots have no de- 
tectable [S 11] emission. Adopting shock velocities of 1530- 
2000 km s" 1 , the distance to Kepler's remnant is estimated to 
be 2.9 + 0.4 kpc. 

4.1.5. RCW86 

Generally — but not universally — accepted as the remnant 
of SN 185, the shell-like X-ray and radio source is found by 
Long & Blair 0] to contain BDSs (W = 660 + 90 km s"\ 
4/4 = 1.0 + 0.2). Figure 3 of their study is particularly illu- 
minating in that it shows sample spectra of radiative and non- 
radiative shocks on the same plot. 

Smith J95ll catalogues a network of B aimer-dominated fila- 
ments surrounding almost the entire remnant and argues that 
they provide a more complete picture of it than the radiative 
filaments do. The morphology of the SNR appears to agree in 
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Figure 4: Combined optical and X-ray image of northeastern RCW 86 (a Galac- 
tic supernova remnant), adapted from Helder et al. 1 39] . The blue and pink 
colors represent the Chandra broad-band (0.1-10 keV) and VLT narrow-band 
Ha emission, respectively. Courtesy of Eveline Helder. (North is up and east is 
left.) 

the Ha, X-ray (from ROSAT) and radio (at 843 MHz) images, 
implying that the non-radiative, optical filaments are a reliable 
tracer of the blast wave. 

Ghavamian et al. [31] measure W = 562 ±18 km s" 1 and 
hlh = 1.18 ± 0.03 from the southwestern part of the rem- 
nant, flanking bright, radiative shocks. Ghavamian, Laming 
& Rakowski 113411 analyze previously unpublished data to add 
two more data points from the eastern (W = 640 ±35 km 
s"\ I b l I n = 1.0 ± 0.2) and northern (W = 325 ± 10 km s"\ 
hlh = 1.06 ±0.10) rims. 

Helder et al. [ 39] use Chandra X-ray proper motion mea- 
surements to measure a shock velocity of v s = 6000 ± 2800 
km s" 1 for northeastern RCW 86. By contrast, the broad Ha 
line has W = 1100 ± 63 km s" 1 . If one assumes « W, then 
the post-shock temperature is T s « 2 keV, in sharp contrast to 
the temperatures of 42-70 keV derived from the proper mo- 
tion measurements, suggesting that a substantial fraction of the 
shock energy is being channelled into cosmic rays (see 38.31) . 

4.2. SNRs in the Large Magellanic Cloud 
4.2.1. SNR 0519-69.0 

Tuohy et al. 1110011 are the first to discover BDSs from SNRs 
outside of our Galaxy. SNR 0519-69.0, a prominent extended 
source in a Ha image but invisible in [O m] imaging, is the only 
remnant in their study for which a reliable detection of the broad 
Ha line is reported, with W = 2800 ± 300 km s" 1 and hlh = 
0.4-0.8. Smith et al. H later report W = 1300 ± 200 km s" 1 
from a slower BDS in the same remnant; they also estimate the 
age of the remnant to be about 500-1500 yrs. 

From the intensity of Ha and Hj3 emission, the density of 
ambient, neutral hydrogen is estimated by Tuohy et al. to be < 



0.1 cm -3 , in contrast to the total ISM density of about 4.7 cm" 3 
based on analyzing X-ray spectra obtained from the Einstein 
Observatory. This suggests that the ISM around SNR 0519- 
69.0 is at least 95% ionized, but such an interpretation rests 
on the assumptions that the X-ray-emitting plasma has cosmic 
abundances and is in collisional ionization equilibrium. If the 
plasma is enriched by the SN ejecta, then the total ISM density 
is instead > 0.3 cm -3 . 

Tuohy et al. also note that SNR 0519-69.0 appears to be a 
weak radio source at 408 MHz compared to other SNRs in the 
Large Magellanic Cloud (LMC), implying that either the ISM 
around it is different or that there exists a fundamental differ- 
ence in the nature of the particle acceleration process in SNRs 
with BDSs. They conclude by speculating that based on the 
similarity of the observed optical, X-ray and radio properties 
with the Galactic remnants of Tycho, SN 1006 and Kepler, the 
progenitor of SNR 0519-69.0 is probably a Type la SN. 

Ghavamian et al. ll35ll use the Far Ultraviolet Spectroscopic 
Explorer satellite to measure broad Lyj3 (W = 3130 ± 155 km 
s" 1 ) and O vi ,11032 (W = 4975 ± 1830 km s" 1 ) line widths; 
the detection of bright (~ 10-20% of Ly/3) O vi emission con- 
firms the prediction of Raymond et al. 117411 . The aperture used 
is comparable to the projected size of the SNR, implying that 
the broad line widths measured have both intrinsic and and bulk 
velocity broadening components. However, the broad H a line 
width of 2800 ± 300 km s" 1 measured by Tuohy et al. ifTooll 
suggests that the Ly/3 width is dominated by thermal broaden- 
ing. In addition to the presence of strong interstellar absorption 
near the center of the Ly/3 line, the O vi and Lyfi lines are heavily 
blended, further adding to the uncertainties in measuring their 
widths. 



4. 2. 2. SNR 0505-67. 9 (DEM L71 ) 

SN R 0505-67.9 is discovered in the study of Tuohy et al. 
1 100]. Based on its somewhat large size (~ 20 pc, assuming 
a distance to the LMC of 55 kpc) and the presence of isolated 
patches of [O m] detected in imaging, they conclude that it is 
older than SNRs 0519-69.0 and 0509-67.5. Like in the case 
of SNR 0519-69.0, its ambient, neutral hydrogen density is es- 
timated to be < 0.1 cm -3 in contrast to the total ISM density 
of 1 cm -3 (based again on assuming cosmic abundances when 
analyzing X-ray spectra). It is not detected in the radio at 408 
MHz. 

Smith et al. ll93ll report the first detection of the broad Ha line 
width (W - 580 ± 70 km s" 1 ), but record only an upper limit 
for the line ratio {hlh ^ 0.7) as they determine the narrow Ha 
line component to be contaminated by radiative shocks. Inter- 
preting the shock velocity to be 300-800 km s" 1 , they use the 
Sedov approximation to estimate the age of the remnant to be 
2^snr/5vs ~ 10 4 yrs, suggesting that the BDSs are transition- 
ing into radiative shocks. Smith, Raymond & Laming I194I1 fail 
to detect the broad Ha emission, which they attribute to poor 
signal-to-noise in their spectra. 

Broad Ly/3 (W = 1135 ± 30, 1365 ± 75 km s" 1 ) and O vi 
,11032 (W = 740 ± 45, 935 ± 125 km s" 1 ) lines are detected at 
two locations by Ghavamian et al. 113 511 . 
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4.2.3. SNR 0509-67.5 

SN R 0509-67.5 is discovered in the study of Tuohy et al. 
iflQOh . Along with SNR 0519-69.0, it is the only other remnant 
to be detected at 408 MHz and also appears to be a weak radio 
source at this wavelength relative to the other LMC remnants. 
This is the only remnant associated with BDSs for which a mea- 
surement of the broad Ha line width has proved to be elusive. 

However, broad Lyfi emission is detected by Ghavamian et al. 
iHl — they measure W = 3710 ± 400 km s" 1 . Unfortunately, 
the FUSE aperture used is essentially the same projected size 
as the SNR. The absence of the broad Ha line precludes a con- 
sistency check for the degree to which the measured broad Ly/3 
line width is contaminated by bulk velocity broadening, imply- 
ing that the measured W value for Ly/3 should be regarded as an 
upper limit. 



4.2.4. SNR 0548-70.4 

SNR 0548-70.4 is discovered in the study of Tuohy et al. 
JlOOll ; imaging reveals patches of isolated [O m] and the rem- 
nant is undetected in the radio at 408 MHz. The first detection 
of the broad Ha line is reported by Smith et al. (W = 760 ± 140 
km s _1 , I b /I n = 1 . 1 ± 0.2) Q ; they also estimate the age of the 
remnant to be ~ 10 4 yrs. 

4.2.5. SNR1987A 

The reverse shock in SNR 1987 A may be the highest- 
velocity BDS ever observed (v, ~ 8000-9000 km s" 1 ) HH, 
but such an interpretation — at least for the "broad" Ha line — 
is controversial. The neutral hydrogen atoms that produce the 
BDS emission are not from the ISM, but are instead from the 
SN ejecta cooled by adiabatic expansion — they stream across 
the reverse shock at ~ 12000 km s" 1 in the frame of the ob- 
server. The hydrogen emission with the greatest Doppler shift 
comes from direct excitation of the atoms, is analogous to the 
narrow line emission and traces out the surface of the reverse 
shock — "surface emission". In the frame of the observer, the 
broad neutrals produced have lower bulk velocities, are associ- 
ated with lower Doppler shifts and appear to originate from be- 
neath the reverse shock surface — "interior emission" 
The measured ratio of interior to surface emission is interpreted 
to be 4 //„ ~ 1 . However, the models of Heng & McCray lEHl 
(see 36.lt predict theoretical values of hlh -0.1 for this range 
of shock velocities, implying that the dominant mechanism re- 
sponsible for the interior emission may not be associated with 
a non-radiative shock without a precursor. 



5. Balmer-Dominated Shocks and the Acceleration of Cos- 
mic Rays 

A number of empirical results suggest that the generation of 
cosmic rays is intrinsic to ~ 1000 km s" 1 BDSs. 



5.7. Narrow Ha Lines with Unusually Broad Widths 

Smith, Raymond & Laming [ 94] obtain high resolution spec- 
tra of the four LMC remnants studied by Smith et al. [93]. With 
an instrumental resolution ~ 13 km s" 1 , they measure narrow 
Ha line widths Wo ~ 30-50 km s" 1 from shocks at multiple lo- 
cations in SNRs 0509-67.5, 0519-69.0, 0505-67.9 and 0548- 
70.4. Such large widths are difficult to reconcile with the pres- 
ence of hydrogen emission, which is not expected to exist above 
~ 10 4 K (Wo ~ 20 km s" 1 ; see §|2]). An important clue provided 
by their observations is that Wo does not appear to correlate 
with v s . Furthermore, several of the spectra show evidence for 
non-Gaussian line wings. 

Smith, Raymond & Laming suggest several mechanisms for 
broadening the narrow Ha lin© 

• Lyman line trapping, in which narrow Ly/3 photons are 
scattered and converted into Ha ones, is a plausible ex- 
planation, but it is found to narrow the line instead (by « 2 
km s" 1 ). 

• Collisions between post-shock protons and pre-shock hy- 
drogen atoms, both prior to and during excitation, may 
heat the atoms and broaden the narrow line component. 
However, both processes are only favorable at high shock 
velocities (> 3000 km s" 1 ), implying that they cannot ac- 
count for the Wo > 30 km s" 1 values measured in the older 
SNRs 0505-67.9 and 0548-70.4. 

• The dissociation of pre-shock molecular hydrogen creates 
neutral hydrogen atoms with velocities ~ 30 km s" 1 , but 
for > 1000 km s" 1 non-dissociative ionization by elec- 
trons dominates over all dissociative channels. Also, H+ 
or H2 are directly converted to fast neutrals via dissocia- 
tion and charge transfer with protons, respectively, without 
ever producing slow neutrals. The ~ 30-40 km s" 1 narrow 
line widths recorded in the young SNRs 0509-67.5 and 
0519-69.0 (which have v s > 1000 km s" 1 ) make it un- 
likely that this is the broadening mechanism at work in 
these remnants. 

• *If temperature equilibration between the electrons and 
ions is nearly complete, then the thermal velocity of the 
electrons exceeds the shock velocity (v±,e ~ 23 v s \ see 
equation (37)), allowing for thermal conduction via the 
electrons to heat the pre-shock gas. Unfortunately, the 
conduction layer is typically orders of magnitude thinner 
than the length scale for ionization, implying that the elec- 
trons are unable to effectively heat the atoms. Photoion- 
ization precursors suffer from a similar difficulty. 

• ^Magnetohydrodynamic (MHD) precursors rely on a low 
pre-shock ionization state, such that the Alfven velocity in 
the ionized gas exceeds the shock velocity upstream. They 
generally require ambient magnetic field strengths of a few 
hundred jiG (e.g., Bi SM ~ 250 jjG for SNR 0519-69.0). 
Furthermore, MHD precursors require a sub-Alfvenic but 



9 In addition to being from a core-collapse supernova. 



10 Precursors are marked with a $ (see Qfor a basic description of precur- 
sors). 
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still supersonic shock — the two components will have 
different shock structures in general (see §3.2 of Draine & 
McKee iB). 

* An interesting possibility that has not been fully explored 
(see 36.2.2b is that a substantial fraction of the broad neu- 
trals, which flow away from the shock front at ~ v s /4, may 
have thermal velocities exceeding the flow velocity. They 
are thus capable of imparting energy to the pre- shock gas. 
This process is expected to be sensitively dependent upon 
the shock velocity, which makes it potentially at odd£3 
with the empirical finding that Wo appears to be uncorre- 
cted with v^. 

*An intriguing possibility is that energetic particles scat- 
tering between the shock front and a pre- shock, turbulent 
region (i.e., cosmic rays) will heat the pre-shock gas (see 
36.2. It . In perpendicular shocks — where the magnetic 
field lines are perpendicular to the flow direction — if the 
broadening is associated with the conservation of the adi- 
abatic invariant 



= constant, 



(15) 



SM 



where v± ~ Wo is the velocity component of an electron 
or ion perpendicular to the field line, then a factor ~ 3 in- 
crease in Wo corresponds to about an order of magnitude 
amplification of the ambient magnetic field. In contrast 
to the other suggested mechanisms, the concern here is 
that the heating is too efficient and the precursor thickness 
should be small enough that the pre-shock gas is not sub- 
stantially ionized. The upper limit to the cosmic ray dif- 
fusion coefficient (equation lfT4l ) is orders of magnitude 
below typical values for the undisturbed ISM, but this can 
be reconciled with the fact that the pre-shock region of a 
BDS is partially neutral — ion-neutral collisions are ex- 
pected to heavily damp the Alfven waves ll23ll . diminish 
the effective value of a:cr, and limit the fraction of shock 
energy deposited into cosmic rays. It is not clear if heating 
by a cosmic ray precursor has a strong dependence on the 
shock velocity. 

Hester, Raymond & Blair 114511 measure Wo ~ 33 km s" 1 at 
multiple locations in the Cygnus Loop, which has the lowest 
shock velocities (v^ < 200 km s" 1 ) among SNRs with BDSs. 
Their kinetic models show that radiative transfer effects are un- 
able to account for this somewhat large width, consistent with 
the conclusions of Smith, Raymond & Laming. 

Ghavamian et al. [ 30] record Wo = 44 ± 4 km s" 1 from 
high-resolution spectra of knot g in Tycho's remnant and also 
detect diffuse Ha emission extending > 1 pc ahead of the knot. 
They attribute the diffuse emission to a photoionization precur- 
sor. Since non-radiative shocks lack recombination zones, they 
reason that the main source of photoionization is photons from 
lines in the immediate post-shock region. Models in which He 11 



11 The velocity-dependence of precursor heating by broad neutrals has not 
been clearly demonstrated ( 36. 2. 2b . 



A304 is the primary source of ionizing radiation from behind the 
forward shock match observed [S n]/Hof and [O m]/H/3 — but 
not [N n]/Hof — line ratios fairly well, but fall short of explain- 
ing the observed narrow Ha line width (To ~ 4 x 10 4 K). They 
speculate that a second precursor might be present. Curiously, 
the narrow Ha line shows weak evidence for the presence of a 
second, broader component — fitting a second Gaussian to the 
line yields a broader component with width ~ 150 km s" 1 — 
the origin of which is unexplained. 

Sollerman et al. 1 96] use high resolution spectroscopy to 
measure narrow Ha line widths in RCW 86 and the remnants of 
Kepler and SN 1006. Their Figure 6 demonstrates that the cor- 
relation between Wo and is not straightforward — the narrow 
line width remains somewhat constant for most of their em- 
ployed data points (Wo ~ 30-40 km s" 1 ), while their highest- 
velocity data point (for SN 1006) has W = 21 + 3 km s" 1 . 
Curiously, SNR 0509-67.5 has W = 25-31 km s" 1 but no 
published measurement of the broad Ha line width, although 
Tuohy et al. H 1 QQfl estimate > 3600 km s" 1 for this remnant. 
Thus, there is tentative evidence for a cut-off in the narrow Ha 
line width at high shock velocities, the origin of which is not 
understood. 

5. 2. Unusually Low Broad-to -Narrow Ha Line Intensity Ratios 

DEM L71 — otherwise known as SNR 0505-67.9 — is the 
poster child for demonstrating that BDS models without pre- 
cursors ( 36.11) fail to capture essential physics. Ghavamian et 
al. 133 1 measure 14 values for W and h //„ around the entire pe- 
riphery of the remnant. While the broad Ha FWHM is < 1000 
km s" 1 , the mean value of hlh is about 0.5; in one location, 
hlh is as low as 0.2. Most of the measured hlh values are too 
low to be explained by either the models of Ghavamian et al. or 
by th e models of Heng & McCray [41] and van Adelsberg et al. 
I lQlh (see ED- 

In addition to broadening the narrow line, precursors con- 
tribute to the narrow line intensity, hence lowering the value of 
hi I n , conceivably explaining the discrepancy between models 
and observations and motivating the need for improved models 
that include cosmic ray precursors (see 36.21) . Ghavamian et 
al. l33ll find no evidence for diffuse Ha emission ahead of the 
BDSs, implying that the precursor thickness is comparable to 
or less than the thickness of the observed filaments. 

Rakowski, Ghavamian & Laming lfv3h confirm the low h lh 
values obtained by Ghavamian et al. with higher-resolution 
spectroscopy, reporting measured values of hlh at 16 posi- 
tions, again with a mean value of about 0.5. 

5.3. Acceleration of Cosmic Rays in Partially Neutral Media 

5.3.1. The Remnant of SN 1006 

It is previously known that the remnant of SN 1006 is "bipo- 
lar" or barrel- shaped, being brightest (in the radio) in the north- 
east and southwest. Early X-ray spectra of SN 1006 are entirely 
featureless, hinting at a non-thermal (i.e., synchrotron) origin. 
Coupled with the absence of evidence for a central pulsar, this 
led to speculation that the non-thermal X-rays originate from 
electrons accelerated at the shock front by the first-order Fermi 
mechanism. 
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Figure 5: Combined optical and X-ray image of the remnant of SN 1006, 
adapted from Cassam-Chena'i et al. 1 16]. Shown are the Chandra synchrotron 
(2-4 .5 keV; in cyan) Hoi and CTIO Curtis Schmidt narrow-band Ha (in red) 
ll09ll emission. Courtesy of Gamil Cassam-Chena'i and Jack Hughes. (North is 
up and east is left.) 



In an insightful study, Winkler & Long [ 108] compare Ha, 
X-ray (ROSAT) and radio (1370 MHz) maps of SN 1006. Pre- 
dominantly thermal X-ray emission is found around the entire 
periphery of SN 1006, spatially offset from the Ha filaments 
and presumably tracing out the post-shock cooling region. In 
the northeast and southwest, the X-ray spectrum is better de- 
scribed by a power law; the detailed correspondence between 
the X-ray and radio structures support a common synchrotron 
origin for emission in both limbs. 

If the ambient magnetic fields are amplified to ~ 100 fiG, 
then radio synchrotron emission at ~ GHz frequencies is asso- 
ciated with E e ~ 1 GeV cosmic ray electrons, 



Vradio ~ 0.5 GHz 



VI GeV/ \ 



_#ISM_ 
100 yUG 



(16) 



while X-ray synchrotron emission at 
from ~ 10 TeV cosmic ray electrons, 



keV energies originates 
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Cosmic ray protons at ~ keV energies are expected to produce 
non-radiative optical emission (see 38.51) . A review of particle 
acceleration in SNRs can be found in Reynolds [ 80] . 

Optical, Ha filaments are detected around much of the pe- 
riphery of SN 1006, but are most prominent in the density- 
enhanced northwest region, which has neutral densities ~ 0. 1-1 
cm -3 compared to ~ 0.01-0.1 cm -3 for the rest of the remnant. 
A sharply-defined radio shock is not seen at the northwest limb; 
together with the lack of correlation with the X-ray emission, 
this suggests the synchrotron component to the emission is neg- 
ligible. This picture is corroborated by the fact that the spatial 
variation of the X-ray surface brightness behind the shock front 
is adequately modelled by a cooling thermal plasma. The spa- 
tial offset between the peaks of the X-ray and Ha intensities 
implies a temporal delay of ~ 120 yrs between the ionization of 



neutral species and the production of thermal X-rays; the radio 
emission develops ~ 600 yrs after ionization. 

Cassam-Chena'i et al. Ill6ll re-investigate the optical, X- 
ray and radio properties of SN 1006, focusing on the poorly- 
explored southeast limb where a somewhat fainter — but still 
clearly-defined — Ha filament can be seen. Again, synchrotron 
X-ray and radio emission appear to be suppressed. However, 
they also find that over small ranges near the eastern and south- 
ern azimuths, synchrotron and non-radiative Ha emission ap- 
pear to be co-spatial (see FigureO, making SN 1006 one of the 
rare cases where such a correlation has been observed. 

Using the BDS in Ha and thermal X-rays as respective trac- 
ers of the blast wave (BW) and contact discontinuity (CD) in 
SN 1006, Cassam-Chenai et al. are able to quantify the ratio of 
the BW to CD radii, 7?bw/^cd- They record /?bw/^cd ~ 1 in 
places where synchrotron emission is present; near the south- 
east azimuth, where synchrotron emission is suppressed, they 
measure Rb w /Rcd ~ 1 • 1 • These empirical values are consis- 
tent with the notion that the production of cosmic rays reduces 
the pressure between the BW and CD, and allows them to be 
located closer to each other (as Warren et al. Ill06ll find for 
Tycho's remnant). Also, it is found that the brighter the BW, 
the smaller the separation between the BW and CD. Using cos- 
mic ray-modified hydrodynamic models, Cassam-Chenai* et al. 
are able to account for the azimuthal variation of the X-ray and 
radio brightnesses by allowing the injection rate of particles, 
ambient magnetic field and cosmic ray diffusion coefficient to 
vary azimuthally. A major shortcoming of their model is the in- 
ability to produce 7?bw/^cd < 1.1 values; an ad hoc reduction 
factor of 12% is applied to all of their results on the assump- 
tion that it is caused by an effect unrelated to the cosmic ray 
acceleration process. 

5.3.2. Tycho's Remnant 

The presence of thin, bright, X-ray synchrotron-emitting 
rims — identified in the Chandra 4.1-6.1 keV banc0 by War- 
ren et al. il06l1 — at the BW all around Tycho's remnant 
implies that there exist locations where BDS and non-thermal 
emission are correlated, although these regions have not yet 
been scrutinized to the same level as for the remnant of SN 
1006. Warren et al. conclude that the mere existence of the 
BDSs argues for the efficiency of cosmic ray acceleration to 
vary around the rim of Tycho's SNR. 

6. Theoretical Models 

Theoretical models of BDSs can be broadly categorized into 
those without ( 36. It and with ( 36. 2b precursors. There is a gen- 
eral disconnect between the models in 36. H and 36.21 — the for- 
mer compute detailed, broad hydrogen line profiles and hll n 
values while ignoring physics that broaden the narrow hydro- 
gen line, while the latter isolate precursor physics without uni- 
fying them with the non-precursor models in a self-consistent 
manner. 



12 Warren et al. estimate the thermal contribution to the 4.1-6.1 keV band at 
the rims to be ~ 10%. 
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6.1. Models Without Precursors 

Observations of BDSs from SNRs commonly yield hlh < 
1 , which early models are unable to account for. For example, 
the models of Smith et al. ll93ll predict hlh > 1, due to various 
simplifications in the way they treat the atomic interactions. 

Ghavamian et al. [31] improve upon the model of Cheva- 
lier, Kirshner & Raymond |19] by including an approxi- 
mate treatment of the shock transition zone — the region 
~ 10 15 (1 cm -3 /ft) cm where excitation, ionization and charge 
transfer take place — and a Monte Carlo calculation of Lyman 
line trapping (see l6.3.lt in their models. The latter involves fol- 
lowing the propagation, absorption and conversion of Lyman 
photons into B aimer ones in the pre- shock region. Electron- 
proton temperature equilibration is parametrized by 



To compare models with observations (of the Cygnus Loop, 
Tycho's SNR and RCW 86), Ghavamian et al. compute hlh 
versus < / eq < 1 for different values of the pre- shock ion frac- 
tion /i on . These calculations are then compared to the observed 
ranges of hlh to constrain the range of allowed / eq values. The 
dependence on f lon is found to be somewhat weak, but in the 
case of southwestern RCW 86 does allow f lon > 0.9 to be ruled 
out. In the case of Tycho's knot g, the Ghavamian et al. models 
are unable to produce the observed hlh = 0.67 ±0.1 value for 
Ha; the theoretical predictions are > 1. 

In the Ghavamian et al. model, the effects of multiple charge 
transfers on the broad neutral distribution are not treated. The 
pre- shock temperature To is a free parameter that is set to 5000 
K. This technique of constraining the degree of electron-proton 
equili brati on is similar to that employed later by van Adelsberg 
et al. rfToTh (see equation ED). 

Heng & McCray 114111 revisit the theory of BDSs and present a 
formalism that allows one to compute the broad neutral veloc- 
ity distributions, rate coefficients and probabilities for atomic 
interactions from first principles, with the only input being the 
atomic cross sections. Multiple reactions involving the same 
atom are visualized using a "reaction tree" (their Figure 4), al- 
lowing for multiple charge transfers and excitations to be in- 
cluded in computing the rate coefficients. Consequently, the 
Heng & McCray model predicts theoretical values for W and 
hi I n as a function of v s and fi\ only the j3 = 0.25 and j3 = 1 
models are published. The narrow hydrogen line is approxi- 
mated by a Dirac-delta function in velocity space, i.e., To - 0. 
The velocity offset between the centers of the narrow and broad 
distribution profiles yields the inclination of the plane-parallel 
shock relative to the line of sight, e.g., « 6° for the Kirsh- 
ner, Winkler & Chevalier [52] observation of the nearly edge-on 
BDS in Tycho's remnant (see 34.1.31) . 

Heng et al. J42I1 include the effects of f{ on and To into the 
Heng & McCray model by treating the shock transition zone in 
the approximation that the broad neutrals and post-shock ions 
share a common bulk velocity and temperature. The expecta- 
tion that the Rankine-Hugoniot jump conditions are obeyed at 
the shock front comes naturally out of the models. The models 
also allow the spatial emissivities of the narrow and broad Ha 



lines to be calculated, providing a consistency checlQ on the 
Raymond et al. 117 811 computations for SN 1006. Lyman line 
trapping is not treated in the Heng et al. model. 

Van Adelsberg et al. II 1 lh generalize the Heng et al. model 
to fully incorporate the effects of multiple charge transfers 
on the broad neutral population and relax the assumption that 
broad neutrals and post-shock ions share a common bulk ve- 
locity. Lyman line trapping is treated using an approximate fit- 
ting function. The more rigorous treatment of the shock tran- 
sition zone produces lower theoretical values of 4/4, thus al- 
lowing self-consistent values of and J3 to be derived from 
measured values of W and hlh (see 38.11) . For example, for 
the W = 1765 ± 110 km s"\ I b /I n = 0.67 ± 0.1 values observed 
in Tycho's knot g ll3lll . van Adelsberg et al. obtain v s « 1600 
kms" 1 and/?;* 0.05. 

6.2. Models With Precursors 
6.2.1. Cosmic Ray Precursor 

Rakowski, Laming & Ghavamian [72] study electron heating 
by lower hybrid waves, which are electrostatic ion waves prop- 
agating nearly perpendicular to the magnetic field lines. They 
have a frequency equal to the geometric mean of the electron 
and ion gyrofrequencies. The growth of lower hybrid waves is 
most efficient in quasi-perpendicular shocks, i.e., the magnetic 
field is nearly perpendicular to the flow across the plane-parallel 
shock. Lower hybrid waves have the property that their phase 
velocity parallel to the field lines greatly exceeds its perpen- 
dicular counterpart, thereby allowing collisionless energy ex- 
change to occur between ions moving across — and electrons 
moving along — the lines. 

Above the critical Alfvenic Mach number of Ma&h ~ 12- 
60, cosmic ray ions channel their energy into amplifying the 
magnetic field; even an initially quasi-parallel shock evolves 
into a quasi-perpendicular one via the development of a helical 
magnetic field. Below Al4,crib the dominant mechanism for the 
cosmic ray ions to deposit energy is through the generation of 
lower hybrid waves, which accelerate a small fraction of elec- 
trons that are in resonance. These accelerated electrons com- 
municate their energy to the rest of the thermal population via 
Coulomb interactions. Magnetic field amplification and lower 
hybrid waves work together — as the field close to a high Mach 
number-shock front gets amplified, the effective value of Ma 
is reduced, eventually allowing lower hybrid wave growth to 
take over. Heating is accomplished via damping of the waves. 
Consequently, there are two precursors: an extended (~ 10 18 
cm) precursor in which amplification occurs, followed by a nar- 
rower (~ 10 15 cm) precursor in which electron heating happens. 

A key assumption of the Rakowski, Laming & Ghavamian 
model is that cosmic ray ions are ubiquitous in SNR shocks. 
If their hypothesis regarding lower hybrid waves is correct, it 
implies that cosmic rays are not by-products of collisionless 
shocks — rather, they are an intrinsic component that affects 
the shock structure, dynamics and energetics. 



3 It is worth noting that for the purpose of predicting the spatial shift between 
the narrow and broad line emissivities — at least for SN 1006 — it is sufficient 
to use the two-component model of Heng et al., which is analytical. 
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Wagner et al. 1 105] use a two-fluid approximation to model 
the back-reaction of the cosmic rays on the thermal, pre-shock 
gas — the cosmic rays are treated as a massless fluid exert- 
ing a bulk pressure; their particle distribution is not followed 
in the calculation. Heating is accomplished by the damping of 
sound waves emitted by the cosmic rays via an acoustic insta- 
bility. Shock parameters appropriate to knot g in Tycho's SNR 
are adopted: v s = 2000 km s _1 , T = 1.2 x 10 4 K, n = 1 cm" 3 
and kqr = 2x 10 24 cm 2 s" 1 . Cosmic ray injection at the gas sub- 
shoclo is included in the model and shown to be necessary in 
order to match the observed spatial emissivities of narrow and 
broad Ha emission. However, matching the spatial emissivities 
also requires a pre-sub-shock temperature of 10 5 K, somewhat 
higher than the To ~ 4 x 10 4 K value implied by the Lee et 
al. ll59ll observations of Tycho's knot g if interpreted within the 
context of a cosmic ray precursor. 

Initially, the shock is weakly modified by cosmic rays, 
but rapidly (~ 1500 yrs) evolves into a cosmic ray-dominated 
shock. A transient solution for the shock that satisfies the ob- 
servational constraints develops within 420 yrs after the initial 
state. About 1000 yrs post-initialization, a density spike devel- 
ops in the solution due to a temporary over-compression associ- 
ated with both the gas sub-shock and the cosmic ray precursor. 
Wagner et al. predict that the density spike may manifest itself 
as a factor ~ 2-4 enhancement in the thermal X-ray intensity. 
Furthermore, they discover a branch of steady solutions with 
low cosmic ray acceleration efficiency — which they do not 
match to the observed spatial emissivities — that are relatively 
insensitive to the shock parameters. Wagner et al. speculate that 
these solutions may account for the lack of correlation between 
Wo and v s (see 35.11) . 

6. 2. 2. Broad Neutral Precursor 

Lim & Raga [60] investigate the broadening of the narrow 
hydrogen line via elastic collisions with broad neutrals. Their 
study considers a range of pre-shock ionization fractions, since 
it is likely the main parameter that controls the upstream energy 
flux, but only for face-on shocks with = 225 km s" 1 . How- 
ever, they find a nearly negligible dependence: for f on = 0.5- 
0.99, Wo decreases from about 22.6 to 21.6 km s" 1 . 

The atomic cross sections and rate coefficients employed in 
the calculations of Lim & Raga are from references in or before 
1970 (see 36.31) . The key process studied, ion-neutral elastic 
scattering, is modelled using an approximate rate coefficient. 
A Gaussian distribution is assumed for the broad neutrals — 
while this is a good approximation at v s = 225 km s" 1 , it breaks 
down at high shock velocities because charge transfer becomes 
inefficient compared to ionization (see 36.31) . It remains to be 
shown how sensitively the increase in Wo depends on v^. 



6.3. The Relevance of Atomic Physics 

6.3.1. Charge Transfer, Excitation and Ionization 

The collisionless, non-radiative shock model [Oil [l9l 41] 
predicts that cold, pre-shock neutral atoms stream past the 



shock front without acceleration or ionization, and subse- 
quently engage in charge transfer/exchange, excitation or ion- 
ization reactions with post-shock electrons and ions. The cross 
section for ionization exceeds that for charge transfer to the 
ground state at just under Av = 3000 km s" 1 (e.g., see Figure 
1 of 114 111 ), somewhat coincident with the thermal velocity of an 
electron with 13.6 eV of energy, v± e ~ 2700 km s" 1 (equation 

ax 

The dominance of charge transfer/exchange for Av < 3000 
km s" 1 (v^ < 4000 km s" 1 ) implies that the broad Ha line pro- 
file and projected^ post-shock proton distribution are virtually 
identical. In this regime, we regard charge transfer as being "ef- 
ficient". The mapping becomes more complicated at higher ve- 
locities because excitation and ionization competes effectively 
with charge transfer. 

Direct collisional excitation of the pre-shock atoms, 



H + {*",/>} ->H* + {<?",/>}, 



(19) 



results in narrow hydrogen emission lines with widths that re- 
flect the thermal conditions of the pre-shock gas (-10 km s" 1 ; 
but see ®. For < 2000 km s" 1 , the rate coefficients are 
dominated by electron impact excitation; for > 2000 km s" 1 , 
excitation by protons takes over (see Figure 3 of 114 ill ). Singly- 
and doubly-ionized helium ions may also excite the atoms, but 
this is a secondary effect unless the abundance of helium is ap- 
preciably more than 10%. 

A secondary population of hydrogen atoms, i.e., the broad 
neutrals, is created when the pre-shock hydrogen atoms ex- 
change electrons with post-shock protons, 



p + H 

/? + H* 



(20) 



where the broad neutrals are labelled as H. Charge transfers to 
both the ground and excited states have to be accounted for in 
any reasonable model 

Experimentally, one can distinguish between elastic scatter- 
ing and charge transfer at high energies because the fast emerg- 
ing particle is either a proton or a hydrogen atom [86]. How- 
ever, quantum indistinguishability prevents such a distinction at 
very low velocities 1 eV); the quantum spin is then the only 
way to label the proton. Consequently, charge transfer at very 
low energies is termed "spin exchange". 

It is best to think ofBDS emission as a probabilistic event — a 
given hydrogen atom may engage in charge transfer, excitation 
or ionization at any moment [41 ]. Its ultimate fate is always 
ionization; the question is how many Balmer, Lyman, Paschen, 
etc., photons will it emit before perishing? 

The BDSs that have been observed so far occur in a velocity 
range that is interesting for atomic physics. The relative veloc- 
ity between interacting atoms and ions, Av ~ 3v^/4 ~ 1000 
km s" 1 , is in the range where neither the processes of charge 
transfer, excitation nor ionization can be ignored. 



14 In cosmic ray physics, the fluid layer where plasma heating occurs is 
termed the "gas sub-shock". 

15 Meaning that the adiabatic index of the gas is y « 5/3 instead of 4/3. 



16 Meaning for a given viewing angle between the shock normal and the line 
of sight. 
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If the pre- shock gas is optically thin to Lyman lines, it has 
Case A conditions. If it is optically thick, then Case B condi- 
tions are present. For example, for Case A about 88% of excita- 
tions to the 3p level result in Lyfi photons (3p — > Is), while the 
remaining 12% yield Ha photons (3p — > 2s); if Case B condi- 
tions prevail, all of the Lyj3 photons are converted into Ha ones 
via scattering, a phenomenon known as "Lyman line trapping". 
Often, the intermediate regime occurs in BDSs, because in the 
case of the narrow Ly (3 and Lyy lines the optical depths are 111911 : 



r(Lyj8) 



r(Lyy) : 



4^(1 -/ion) / Tp r 1/2 

/ion U0 4 K/ 
1.4^(1 -/ion) / To 



/io 



\10 4 K/ 



(21) 



where / on is the pre- shock ion fraction and S is fraction of neu- 
tral gas at a temperature To. The larger number of scatterings 
experienced by Ly/3 photons results in a more complete conver- 
sion of Ly J3 to Ha than of Lyy to HJ3, thus enhancing the relative 
intensity of Ha. As such, the "Balmer decrement", defined as 
the intensity ratio of Ha/Hj3, generally increases with optical 
depth. High values of the Balmer decrement may indicate the 
increased importance of collisions, but this situation is degen- 
erate with the presence of dust, which preferentially attenuates 
HJ3 photons. 

6.3.2. References for Atomic Cross Sections 

For convenience, we summarize the most modern list of ref- 
erences by atomic process: 

• Charge transfer/exchange: Janev & Smith [48]; Belkic et 
al. fi; Harel, Jouin & Pons JH, Schultz et al. & 

• Excitation: Janev & Smith H; Bubelev et al. Ell ; 
Balanga, Lin & Feautrier 



• Ionization: Janev & Smith [48]. 

For processes involving neutral helium, singly-ionized he- 
lium and alpha particles: 

• Charge transfer/exchange: Barnett et al. [0]; Janev & 
Smith H; Hose H. 

• Excitation: Janev & Smith [0]; Bray et al. ifioll. 



|; Angel et al. Ol; 
; Rinn et al. llH ; 



• Ionization: Peart, Walton & Polder 1168 
Rudd et al. H; Shah & Gilbody JM 
Shah et al. EM Janev & Smith H. 



Other references may be found in Laming et al. [58]. Several 
fitting functions for some of the atomic data (in the references 
listed above) are given in the appendices of Heng & Sunyaev 
Hi and van Adelsberg et al. lIlOlll . We note that Schultz et 
al. computed cross sections for proton-hydrogen atom charge 
transfer (and spin exchange) for 200 < v s < 20000 km s" 1 
using a variety of techniques. 



6.3.3. The Scarcity of Atomic Cross Sections 

Even with modern techniques, the uncertainties in the com- 
puted cross sections can be ~ 10-30%, representing a non- 
negligible source of systematic error for any model of BDSs 
that incorporates the effects of atomic physics. In particular, 
the proton-hydrogen atom collision system remains an active 
area of research J28, 62, 86|,|24]]. Cross sections for the excita- 
tion of hydrogen atoms by protons at the velocities relevant for 
BDSs remain scarce. 

The proton-hydrogen atom collision system is qualitatively 
described by three different energy regimes with no clear 
boundaries between them B28I1 . In the low-energy regime, 
charge transfer is the dominant, resonant channel, since al- 
most no energy needs to be transferred from the relative nu- 
clear to the electronic motion — the system may be regarded as 
a quasi-molecule with a temporally-varying internuclear sepa- 
ration [ 86]. In the high-energy regime, charge transfer becomes 
a very weak process compared to excitation and ionization — 
the projectile then provides a small perturbation to the wave- 
function of the target. In the intermediate-energy regime, the 
interplay between charge transfer, excitation and ionization is 
important; in this regime, obtaining detailed results like the 
population of projectile or target sub- shells is computationally 
challenging. 

Roughly speaking, for projectiles and targets of nearly equal 
mass, atomic cross sections are somewhat difficult to compute 
when 

Av m D Av 2 

„ (22) 



1 



a FS c 



25 keV, 



where a^s ~ 1/137 is the fine structure constant. For exam- 
ple, the one-center, close-coupling expansion ll62ll commonly 
used to solve the time-dependent Schrodinger equation, in or- 
der to obtain proton-hydrogen atom excitation cross sections, is 
inaccurate below impact energies ~ 20-30 keV, because charge 
transfer becomes the dominant channel and computational con- 
vergence is thus difficult to achieve. Above energies of a few 
hundred keV, perturbative methods may be employed W28I1 . 

A crude scaling law to obtain excitation — or charge transfer 
to excited states — cross sections to the level n^h is 



0~n 2 l 2 - I — I °~nih, 

w 



(23) 



where o~ nx i x is the cross section to the atomic level n\l\ (n\ < 
ni). This relation should only be used for order-of-magnitude 
computations or as a last resort in the absence of available 
atomic data. 

The glaring lack of cross sections for the excitation of hydro- 
gen atoms by protons to the levels 4s, 4p, 4d, 4/ and above, 
at Av ~ 1000 km s" 1 , implies that Balmer decrements cannot 
be computed accurately enough for them to be used as shock 
diagnostics at v s ~ 1000 km s~ l . 



17 The adjective "coupling" comes from having to solve a pair of coupled 
Schrodinger equations. 
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Figure 6: Montage of narrow-band Ha images of the optical bow shocks asso- 
ciated with pulsar wind nebulae. On disnlay are J2124-3358 (H, B0740-28 
lU, B1957 +20 (H and B2224 +65 Gj]. The arrows mark the positions of 
the pulsars. Courtesy of Bryan Gaensler. (North is up and east is left.) 



7. Balmer-Dominated Shocks in Other Astrophysical Ob- 
jects 

The abundance of hydrogen present in the universe and the 
existence of BDSs across a wide range of velocities suggest that 
they should be present in other astrophysical objects as well. In 
this section, I comment on the implications of both probable 
and possible detections of BDSs. A basic but key test is to 
demonstrate that the detected hydrogen line can account for the 
measured luminosity assuming plausible values for the physical 
parameters, despite originating from a non-radiative shock. 

7. 1 . Pulsar Wind Nebulae 

Old (> 10 9 yrs), milli- second pulsars — neutron stars that 
are spun up by the accretion of material from a companion — 
lose their rotational energy by driving a wind into the ISM. One 
naturally expects the wind to manifest itself as a BDS in Ha, 
but such detections are rar^H because the ISM around these 
pulsar wind nebulae is expected to be either fully-ionized or of 
low density. Nevertheless, several examples have been found 
(Figure [6]), spanning a range of distances (~ 0.1-2 kpc) and 
spin-down luminosities (~ 10 33 -10 35 erg s" 1 ). 

Kulkarni & Hester [ 56] discover, via Ha imaging, a distinct, 
comet- shaped nebula around the pulsar PSR B 1957 +20, which 
has a period of 1 .6 ms. The pulsar appears to be moving through 
the ISM at ~ 100 km s" 1 , producing a shock of comparable 
velocity. About 1-10% of the rotational luminosity (of about 
3 x 10 35 erg s" 1 ) is channelled into Ha emission. The spectrum 
reveals weak or absent [O m] and [S n] emission alongside a 



prominent, narrow Ha line, consistent with the properties of a 
BDS. The broad Ha line is not detected. The narrow H/3 line is 
barely seen — assuming a positive detection, the dust-corrected 
B aimer decrement is > 6, at least twice as large as the expected 
value of 3. 

With improved spectroscopy, Aldcroft, Romani & Cordes ^\ 
clearly detect Ha and Hj3 lines from the bow shock around PSR 
B 1957 +20; the nebular lines of heavy elements have upper lim- 
its of 3.5% of the Ha line intensity. The mean B aimer decre- 
ment is determined to be Ha/Hj3 = 3.5 + 0.5, implying that 
the visual extinction is Ay < 1. A single Gaussian fit to the 
Ha line indicates that the face-on, bulk velocity of the flow is 
-120 km s" 1 . A five-parameter kinetic model is fitted to the 
observed spectrum: the ambient density is 0.41+0.02 cm -3 and 
the distance to the pulsar is 1.2 + 0.01 kpc. Aldcroft, Romani & 
Cordes remark that, at least in the case of PSR B1957 +20, the 
analysis of the Ha line provides the cleanest measurement of 
the distance, as noise in the radio pulsar and its faintness render 
conventional parallex methods unlikely to succeed. 

Following the discovery of the optical nebula around PSR 
B1957 +20, several other examples are found. Cordes, Ro- 
mani & Lundgren 112 ill present a B aimer-dominated spectrum 
of the Guitar nebula produced by the high- velocity (> 800 km 
s" 1 ) pulsar PSR B2224 +65. Bell et al. |8] report an image 
of the Ha bow shock around PSR J0437-4715. Van Kerkwijk 
& Kulkarni II 10311 present spectra and Ha images of the neb- 
ula around the neutron star RX Jl 856.5-3754. Jones, Stappers 
& Gaensler [49] image an optical bow shock associated with 
PSR B0740-28. Pellizzoni, Mereghetti & De Luca image 
a Ha nebula possibly associated with the central X-ray source 
in SNR G266.2-1.2. Gaensler, Jones & Stappers 112911 image 
the Ha bow shock around PSR J2 124-3 35 8. It is worth noting 
that the shocks detected in these pulsar wind nebulae comprise 
both non-radiative and radiative components, based on the ab- 
sence/presence of lowly-ionized metal lines. 

The transverse velocity of the pulsar contributes to its ob- 
served period derivative and affects its derived age and mag- 
netic field strength [8]. The separation between the apex of 
the bow shock and the pulsar is known as the "stand-off dis- 
tance", denoted by r w . It arises from a balance between the ram 
pressure of the ISM (~ ^H^Vp ulsar ) and the pressure from the 
relativistic pulsar wind (~ E/47rcrfy, where v pu i sa r is the pulsar 
velocity and E is its spin-down luminosity. By measuring v pu i sa r 
and r w , the ambient density can be estimated, 



(24) 



4nmncrlvf 



pulsar 



18 Cordes, Romani & Lundgren |2l]] estimate that 
may be able to excite Ha nebulae. 



5000 Galactic pulsars 



For example, Jones, Stappers & Gaensler estimate n ~ 0.25 
cm" 3 forB0740-28. 

Despite establishing the presence of BDS emission in these 
optical bow shocks, the broad Ha line has never been robustly 
detected in any pulsar wind nebula. 

7.2. High-Redshift Galaxies 

Heng & Sunyaev 0] reason that the same ~ 1000-10000 
km s" 1 BDSs — produced by SNRs in the Sedov-Taylor phase 
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— that produce broad B aimer lines should also generate lumi- 
nous (by factors ~ 10-100) Lya lines in young, high-redshift 
galaxies. The broad wing^J of the Lya line allow it to escape 
the local optical depth of the SNR, as demonstrated by the Ly/3 
detections of BDSs by Ghavamian et al. 113511 . At redshifts of 
z ~ 2, the ultraviolet Lya line is red-shifted into the optical, thus 
allowing it to avoid Galactic extinction. If multiple B aimer and 
Lyman lines are detected from these distant galaxies, they can 
be used to estimate the average shock velocities and J3- values of 
the constituent SNRs. 

Shapiro et al. [92] use high-resolution integral field unit 
spectrographs on 8-10 m-class telescopes to obtain near- 
infrared spectra of 47 galaxies at redshifts of z ~ 2. Four of 
these galaxies have previously-identified active galactic nuclei 
(AGN). Stacking the spectra of the (putative) non-AGN galax- 
ies reveals a broad component to the Ha line with W > 1500 
km s" 1 ; in comparison, the 4 galaxies with AGN have a broader 
(~ 3000 km s" 1 ) and stronger (by a factor ~ 3) Ha line. In 
stacked spectra from 7 well-resolved galaxies, the strength of 
the broad Ha line is observed to be greater near the galactic 
centers (< 3 kpc) than in the outer regions (3-15 kpc). 

By grouping the 47 galaxies into three stellar mass bins 
(Af* < 2xlO lo M , M* = 2-7xlO lo M andM* > 7xlO lo M ), 
Shapiro et al. demonstrate that the strength and width of the 
broad Ha line is correlated with M*, which can be interpreted 
in two ways. Higher stellar masses stem from higher star forma- 
tion rates — it is plausible that galactic winds driven by ~ 10 4 
SNRs in each galaxy are powering the broad Ha emission via 
the creation of broad neutrals in BDSs, as suggested by Heng 
& Sunyaev. Such a mechanism is energetically consistent with 
the measured luminosities (~ 10 41 -10 42 erg s" 1 ) because the 
ISM in these z ~ 2 galaxies is ~ 10-30 times denser, thereby 
boosting the broad Ha luminosity by factors ~ 100-1000. 

The second possibility is that the broad Ha emission orig- 
inates from the broad line regions around unobscured AGN. 
Higher stellar masses imply higher bulge and (possibly) black 
hole masses. By using an empirical mass-luminosity relation to 
estimate the black hole masses (to within a factor ~ 3), Shapiro 
et al. show that the black holes are about an order of magnitude 
less massive (relative to their stellar bulge masses) than local el- 
liptical and bulge-dominated spiral galaxies. This is interpreted 
as evidence for rapid bulge formation in z ~ 2 galaxies preced- 
ing a few Gyr of super-massive black hole assembly. 

The data obtained by Shapiro et al. do not allow one to ro- 
bustly distinguish between these two possibilities. 

7.3. Symbiotic Stellar Systems 

Symbiotic stars are binary systems typically comprising a 
white dwarf (WD) and a cool giant star. Both stars lose mass via 
winds and the colliding winds lead to complex spectral behavior 
at all wavelengths from the radio to X-ray. CH Cygni is perhaps 
the most studied of symbiotic stars [ 20] . In 1977, it produced a 



19 Heng & Sunyaev termed these contributions to be "non-thermal", but this 
is not meant to be associated with particles that produce non-thermal X-ray 
and radio emission. Rather, it is to indicate hydrogen line emission that is not 
narrow and typically associated with recombination. 



powerful outburst that lasted until 1986 — bipolar radio and op- 
tical jets appeared towards the end. A mysterious Lya line with 
width > 4000 km s" 1 appeared in 1985; during the quiescent 
phases that followed (1990-1991, 1995), it disappeared along 
with the hot ultraviolet and optical continuum, then re-appeared 
during the 1992-1995 active phase (albeit with its line strength 
weakened by a factor ~ 2-3). Ultraviolet spectra from the In- 
ternational Ultraviolet Explorer show that the Lya line is by 
far the most prominent, outshining the singly- to triply-ionized 
metal lines by at least an order of magnitude. 

Contini, Angeloni & Rafanelli [20] suggest that broad Lya 
emission is generated in a BDS from the WD that propagates 
away from CH Cygni (as opposed to the WD wind that collides 
with the wind from the giant star). The contemporaneous ap- 
pearance of the [O m] /15007 line limits the ambient density to 
n < 10 4 cm -3 in order to be consistent with avoiding collisional 
de-excitation. If the [O m] and Lya emission have the same ori- 
gin, then they cannot be from the region in between the WD 
and giant star, which has densities ~ 10 8 -10 9 cm -3 . 

Selvelli & Hack [ 89] measure the Lya flux to be about 2. 16 x 
10" 10 erg cm -2 s" 1 . Since the distance to CH Cygni is 270 pc, 
the Lya luminosity is about 2 x 10 33 erg s" 1 . The theoretically- 
expected value is 

« 47rR BW n (1 - /ion) v s e bL y a E Lya , (25) 

where Rbw ~ 2.25 x 10 17 cm is the blast wave radius and 
n{\ - /ion) ~ 0.2 cm -3 is the pre-shock neutral density lEoll . 
The energy of a Lya photon is E^ ya « 10.2 eV, and we take 
the number of broad Lya photons produced per ionization to be 
ebLya ~ 1 . If we take the reported lower limit for the line width 
to be the shock velocity (without correcting for bulk velocity 
broadening), then we get L^ ya « 8 x 10 32 erg s" 1 . Therefore, 
the broad Lya line detected in CH Cygni may be energetically 
consistent with being from a BDS, although further investiga- 
tion is required. 

7.4. Supernovae 

Komossa et al. l54ll detect Ha and Hfi lines with com- 
plex structures from a transient event in the non-active galaxy 
SDSSJ095209.56+214313.3. These line complexes are respec- 
tively fitted by 7 and 5 Gaussians: 1-3 components for the nar- 
row Ha and [N 11] /16548, 6584 lines; 2 narrow "horns"; and 2 
broad Ha line components. For both Ha and HJ3, they report 
W « 2100 km s" 1 ; the line intensities relative to [O m] A5001 
are 9.7 and 0.79, respectively. 

Using the ratio of the [S 11] /16716, 6731 doublet, the ambient 
density is estimated to be n ~ 10 2 cm" 3 . However, the [O m] 
/14363,5007 doublet yields n ~ 10 7 cm -3 , clearly indicating 
that the oxygen lines are originating from a different region. 
Adopting n > 10 2 cm -3 implies that the size of the region as- 
sociated with narrow Hj3 emission is < 50 pc. 

It is uncertain if this transient event is indeed a SN. However, 
one can describe some characteristics of this event if interpreted 
in the context of a BDS associated with a SN explosion. Mys- 
teriously, the Balmer decrement for the broad lines is > 5.7 
and goes up to about 12.2, higher than that for the narrow lines 
(2.6-3.2); the narrow Balmer decrement is consistent with Case 
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B recombination Il67h . In the contexts of both dust extinction 
and Lyman line trapping, it is somewhat difficult to construct 
scenarios where the broad lines are more susceptible than the 
narrow ones. 

If all of the broad Ha emission is due to non-radiative shocks 
of age t agQ , then using equation © to assert t YQC > t agQ implies 



n < 5 x 10 6 cm -3 



10 yrs 



-1 / \0.7 

m, 



m, 



2100 km s 



-l 



, (26) 



if we assume ~ W. If it originates from a BDS, the broad Ha 
is emitted at a distance from the SN of 



^bds ^ 



rec^age 



a 



2 pc, 



(27) 



where Lbu a ~ 2 x 10 41 erg s" 1 is the mean broad Ha luminosity 
observed 115411 . We take the number of broad Ha photons per 
ionization to be €bna ~ 0.2. 

Indeed, Ha lines with widths ~ 1000 km s" 1 are commonly 
detected in supernovae, but the broadening mechanism is usu- 
ally assumed to be bulk fluid motion. A rare (< 1% by detected 
number) and intriguing sub-class are the "Type Ia/IIn" super- 
novae that exhibit spectral features typically associated with 
thermonuclear explosions alongside prominent hydrogen lines. 
A striking example is SN 2005gj, which shows both narrow and 
broad components of Ha, HJ3, Hy and He i /15876, 7065 lines 
H. The broad Ha and Hf3 lines have widths of W ~ 2000 km 
s" 1 and luminosities ~ 10 40 erg s" 1 . 

It is tempting to interpret these broad hydrogen lines in the 
context of BDSs, in part because the lines are observed to be 
symmetric, somewhat Gaussian and free of distortions typi- 
cally associated with geometrical effects (e.g., flat-topped or 
double-horned profiles); the measured B aimer decrements are 
also large (~ 6). However, the required mass of the circumstel- 
lar medium is prohibitive: 



M CS M 



UOMrr 



I 0.2 ) 



(28) 



where Ebna ~ 10 47 erg is the energy of the broad Ha line 
released over ~ 100 days, as observed by Aldering et al. 
Since Balmer emission is recorded from 11 to 133 days post- 
explosion, the volume of the circumstellar material is ~ 10 50 
cm 3 if the ejecta velocity is taken to be 25000 km s" 1 , implying 
that the ambient density is n ~ 10 9 cm -3 . Therefore, even if 
shocks are responsible for producing the broad Ha lines in SN 
2005gj, they are probably radiative. 

7.5. TTauri Stars 

T Tauri stars — young, luminous and pre-main sequence — 
have been observed to emit symmetric Ha lines with broad 
wings often extending to about 500 km s" 1 , somewhat higher 
than expected estimates for the infall velocity of gas (200-300 
km s" 1 ). The Ha lines typically have a blue-shifted absorption 
trough near line center. 

Ardila et al. 0] detect Ha lines with widths of W « 120-330 
km s" 1 in DR Tauri, DF Tauri, BP Tauri and RW Auriga (listed 
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Figure 7: Ratio of electron to proton temperatures immediately behind colli- 
sionless shock fronts as a function of the shock velocity, obtained via analysis 
of the two-component Ha lines from B aimer-dominated shocks around super- 
nova remnants. Based on calculations from van Adelsberg et al. 110111 . 



in order of ascending W). For these objects, Mg n A2196, 2803 
lines are observed to have spectral profiles that are virtually co- 
incident with Ha, suggesting a common origin for these lines 
and ruling out the linear Stark effect as an explanation (which 
predicts narrower Mg n than Ha lines) . Curiously, the HJ3 lines 
— which are only detected in DF Tauri and BP Tauri — have 
widths of about 130-170 km s" 1 , somewhat lower than the Ha 
line widths. No correlation is found between the accretion rate 
and the Mg n line fluxes. Balmer decrements are not reported 
by Ardila et al. The correlation between the Ha and Mg n lines 
led Ardila et al. to suggest a broadening mechanism that is in- 
dependent of atomic physics, such as Alfvenic turbulence. 

The association between the broad Ha lines and BDSs is sug- 
gestive but unestablished. If true, one can study if there exists 
a correlation between the shock velocity and the abundance of 
crystalline grains that may have been processed by shock heat- 
ing. It is unclear if Mg n emission (which is from a resonant 
doublet) can be produced in such a BDS. 

7.6. Extragalactic SNRs 

The optical discovery of individual, extragalactic SNRs with 
BDSs is unlikely, primarily because non-radiative hydrogen 
emission is very faint and difficult to detect. Past searches 
for distant SNRs were based on the selection criterion of large 
[S n]/Hof ratios, which violates one of the stated properties of 
BDSs (see £Q). 

8. The Future 

8.1. Electron-Ion Temperature Equilibration in Collisionless 
Shocks 

Ha observations of BDSs present a unique opportunity to 
quantify the degree of electron-proton temperature equilibra- 
tion in the immediate post-shock regions, potentially yielding 
insight into the plasma physics of these shock fronts. 

Ghavamian, Laming & Rakowski IK34I1 use 10 measurements 
of W and hll n from four SNRs to derive values of J3 and v s 
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using the model of Ghavamian et al. 113 lh (see 96. lb . For their 
fifth SNR, DEM L7 1 , they use the measured W to infer a range 
of shock velocities and T p — J3 is obtained when combined with 
a Chandra X-ray estimate of T e . The curve 



fi- 



ll v s <400kms -1 , 
£v7 2 v, > 400 km s" 1 , 



(29) 



appears to fit the data fairly well; they also assume min{/3} « 
m e /m p . 

The Ghavamian, Laming & Rakowski model (later expanded 
upon by Rakowski, Laming & Ghavamian [72]) proposes using 
lower hybrid waves (see 96.2b to heat the electrons and predicts 
J3 oc v~ 2 , which implies that the electrons are heated to a tem- 
perature of 

T e =p p *—£-*03teV, (30) 
16k 

independent of the shock velocity, suggesting a process occur- 
ring within the pre- shock medium rather than the shock front 
itself. 

Detailed modelling by van Adelsberg et al. lIlOlll shows that 



input: {W,I b /I n }-+ output: {v„£}, 



(31) 



with a somewhat weak dependence on the pre- shock ion frac- 
tion /i on (which they set to 0.5 for their reported results). Us- 
ing this approach, they are able to self-consistently derive 14 
values of ft = fi(v s ) from 7 SNRs (Figure [7]), three more than 
Ghavamian, Laming & Rakowski. Fitting a j3 oc v~ 2 curve to 
the derived values yields a reduced^ 2 value of about 5, albeit 
with ~ 10-30% uncertainties in the atomic cross sections (see 
96.31) . Also, they find min{/3} « 0.03 ~ ^Jm e /m p as opposed to 
m e /m p . Curiously, their derived ft values hints at a rising trend 
for > 1000 km s" 1 , which remains unexplained. 

We note that the low values 0.06-0.07, derived for the 

northwestern filament in SN 1006, are consistent with the find- 
ings of Laming et al. ll58ll (see 94.1.21) . Rakowski, Ghavamian 
& Hughes flVOfl combine the analysis of optical Fabry-Perot 
spectra (obtained by Ghavamian et al. 113 311 ) with Chandra X- 
ray data to arrive at a similar conclusion for DEM L7 1 — their 
planar shock model yields J3 « 0.01-0.46 at 4 locations with 
faster shocks (v s * 800-1000 km s" 1 ). 

Both the studies of Ghavamian, Laming & Rakowski and van 
Adelsberg et al. fail to accomodate the low hlh values mea- 
sured in DEM L7 1 (see 95.2b and are unable to extract j3 values 
from most locations in this remnant. However, they already 
rule out the model of Cargill & Papadopoulos [15], who sim- 
ulate the heating of the pre- shock medium by shock-reflected 
ions; their models predict j3 ~ 0.2 independent of v s . It is fair 
to say that the technique of quantifying electron-proton temper- 
ature equilibration from optical studies of BDS is still nascent 
and requires further refinement, particularly with respect to in- 
corporating percursor physics. 

The broad O vi line widths measured by Ghavamian et al. 
ll35ll allow an estimate of the oxygen ion-proton temperature 
equilibration, since 



(32) 




where W is the width of broad Lyfi in this case; Wo and To 
are the width of the O vi line and the post-shock temperature 
of the O vi ion, respectively. For SNRs 0505-67.9 and 0519- 
69.0, they measure W/Wo < 4, implying incomplete ion-proton 
equilibration. From three estimated values of To/T p , they ten- 
tatively conclude that there is less ion-proton equilibration at 
higher shock velocities. The far-ultraviolet lines alone preclude 
an estimate of the electron-ion equilibration, because interstel- 
lar dust extinction removes the narrow line component. 

Other studies using Far Ultraviolet Spectroscopic Explorer 
detections of the A 1032, 1038 lines have been performed. Ko- 
rreck et al. JH measure W = 2100 ± 200 km s" 1 from 
SN 1006. Combined with the Ghavamian et al. 113 ill value of 
W = 2290 ± 80 km s" 1 for Ha, this implies T /T p « 13, less 
than mass-proportional heating and indicative of poor equilibra- 
tion. By contrast, Raymond et al. 1 77] infer that the protons and 
oxygen ions in northern Cygnus have the same temperatures to 
within a factor of 2.5. 

The microphysics of electron-ion temperature equilibration 
at collisionless shock fronts is relevant to Sunyaev-Zel'dovich 
studies of galaxy clusters. The Compton y-parameter, inte- 
grated over the area of the cluster, is a proxy for the cluster 

mass M c i uster , 



n P T P dV c 



cluster 



!F(Af c i uster ), 



(33) 



where n e is the electron density, V c iuster is the volume occupied 
by the cluster and T - !F(M c i uster ) is a power-law function of 
the cluster mass. Models that compute T e typically assume 
complete equilibration between the electrons and ions. Rudd 
& Nagai 1 84] show that if partial equilibration is assumed, then 
the value of M c i uster inferred is about 5% higher)^ More impor- 
tantly, the uncertainty in the temperature equilibration translates 
to < 30% uncertainties in the electron temperature as a func- 
tion of distance r c i uster from the center of the cluster, introduc- 
ing an additional source of systematic uncertainty for models 
attempting to extract the position of the virial shock front from 

Te(? cluster)- 

Finally, we note that a review of electron-ion equilibration in 
collisionless shocks from SNRs can be found in Rakowski 117 1R. 



8.2. Test Beds for Particle Acceleration in Partially Neutral 
Media 

The general absence of non-thermal emission at the locations 
where BDSs are present is consistent with the picture that ion- 
neutral interactions damp out the Alfven waves necessary for 
particle acceleration 1E3I1 . but the dependence of this process on 
the physical conditions of the shock needs to be further quan- 
tified. In cosmic ray-dominated shocks, it is possible that suf- 
ficient electron heating is generated such that no neutral atoms 
survive to cross the shock, hence rendering ion-neutral damp- 
ing ineffective/irrelevant. BDSs represent a unique opportunity 
to test this hypothesis in more detail. 



20 This uncertainty is comparable to the ~ 7-8% scatter inherent in using 
simulations to calibrate the relation in equation {33V 
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More interestingly, the locations in the remnants of SN 1006 
and Tycho where non-thermal and BDS emission are correlated 
serve as valuable test beds for models of turbulently-amplified 
magnetic fields in partially neutral media HUES], and therefore 
deserve further scrutiny. 

8.3. Measuring Injection Energy from Face -On Balmer- 
Dominated Shocks 

Raymond et al. Ip74ll and Raymond J75ll suggest a search 
for broad Ha emission from face-on BDSs as a way of testing 
models for shock acceleration of cosmic rays. By comparing 
the velocity shift ~ 2v s /(y + 1) between the narrow and broad 
hydrogen lines with W, one is in essence comparing and T t , 
respectively. This provides an estimate of the energy channelled 
into cosmic ray ions. Specifically, the thermal velocity of the 
protons is reduced by a factor 1 / V 1 + Pcr/Pg, where Pcr and 
Pq are the cosmic ray and gas pressures, respectively B79I1 . For 
example, W is expected to be ~ 30% lower if half of the energy 
goes into non-thermal particles. Quantifying the presence of 
cosmic ray ions allows us to test the Rakowski, Ghavamian & 
Laming [ 72] precursor model of lower hybrid wave heating (see 

SE3. 

Such face-on shocks are extremely difficult to observe, be- 
cause their brightness is not enhanced by limb brightening like 
in the case of edge-on shocks. To date, no reliable catalogues 
of face-on BDSs exist. However, recent work by Helder et al. 
l39ll combines X-ray proper motion measurements with a de- 
tection of the broad Ha line from an edge-on BDS in RCW 86 
to infer that a substantial fraction of the shock energy is being 
channelled into cosmic rays. By empirically constraining the 
parameter space of €qr versus wqr — respectively the fraction 
of escaping energy associated with cosmic rays (relative to the 
kinetic energy of the shock) and the ratio of non-thermal to the 
total post- shock pressure — they infer wqr > 50%. 

8.4. B aimer-Dominated Shock Models with Cosmic Ray Pre- 
cursors 

There is a lack of models for cosmic ray-mediated shocks in 
partially neutral gas where / gyro < Iqr < ^atomic and v s ~ 1000 
km s" 1 . The observations of unusually small and large values 
of hi I n and Wo, respectively, in several SNRs — as well as 
the correlation between optical and non-thermal X-ray emission 
observed in SN 1006 — provide incentive for constructing these 
models. 

8.5. Ha Line Profiles as Tracers of the Post-Shock Proton Dis- 
tributions 

Since charge transfer provides a direct mapping between the 
post-shock protons and broad neutrals, high-precision measure- 
ments of the broad Ha line can be used to quantify the post- 
shock proton velocity distribution, particularly if it deviates 
from being in a Maxwellian distribution. 

Raymond, Isenberg & Laming Ir79ll study the effects of "pick- 
up ions" — scrutinized extensively in the solar wind — which 
in the context of BDSs are protons created when neutral atoms 
pass through the collisionless shock front, become ionized 
downstream and suddenly "see" the ambient magnetic field. 



The protons gyrate around the field lines if the associated time 
scale, t gyro ~ 27rm p c/eBisM, is less than the proton-proton tem- 
perature equilibration time t pp , 

B 1SM > 0.33 yuG ( ^ r ) (-A^ \ (34) 

\ 1000 km s" 1 / \ lcm- 3 10/ 

where we note that t pp = ^/m e /m p t ep with t ep being the time 
scale for electron-proton temperature equilibration (equation 

G2). 

The gyration of these protons forms an unstable ring in 
three-dimensional velocity space, which then evolves towards 
isotropy via scattering and self-interaction. In their study, Ray- 
mond, Isenberg & Laming assume that the end result is a 
bi-spherical distribution II 1 0711 . Pick-up ions introduce non- 
Maxwellian components to the proton distribution function — 
specifically, the contributions occur near the broad Ha line cen- 
ter. Assuming that the broad Ha line exactly traces the proton 
distribution — an assumption that breaks down for > 2000 
km s" 1 shocks because charge transfer is no longer the domi- 
nant atomic process (see 36.3b — Raymond, Isenberg & Lam- 
ing find that the shock velocity can be under-estimated by as 
much as 15%. 

If the magnetic field is laminar, then the broad Ha line is 
generally asymmetric about line center and is offset by a veloc- 
ity shin 21 ! from the narrow Ha line center. If 0b is the angle 
between the flow across a plane-parallel shock and the mag- 
netic field lines, then the shift varies from zero at 0b =0 to 
~ va at #b = 90°, where va ~ 10" 3 -10" 2 v^ (equation Oil) is 
the Alfven velocity. If the magnetic field is turbulent, then the 
broad Ha line is a smeared-out composite of asymmetric lines 
with a range of 0b values — one can in principle extract infor- 
mation about the nature of the turbulence from careful analysis 
of the line. Turbulent diffusion dominates thermalization when 
the associated length scale is 

^^Mi^Hi^f- (35) 

Generally, both scenarios imply that information on the mag- 
netic field geometry can be obtained from high-precision stud- 
ies of the broad Ha line. 

The Raymond, Isenberg & Laming study relaxes the com- 
mon assumption that newly-created protons are instantly as- 
similated into the thermal population. This assumption is es- 
pecially suspect when the pre- shock neutral fraction is high and 
the energy of the thermal, post-shock protons is non-negligibly 
reduced by the act of assimilation Il60n . For example, the pre- 
dicted non-Gaussian contributions to the broad Ha line will be 
much easier to see in Tycho's SNR (/i on « 0.15) rather than in 
the remnant of SN 1006 (/i on ~ 0.9). Since the broad Ha line 
is emitted from a region of varying pick-up ion fraction, the 
creation of pick-up ions should be included in future models of 
BDSs. 



21 This has the implication that edge-on shocks may have a finite velocity 
shift between the centroids of the narrow and broad Ha lines, and is degenerate 
with the interpretation that the shift yields the angle between the line of sight 
and the vector normal to the plane of the shock. 
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If cosmic ray protons are present in the distribution, it is con- 
ceivable that they will introduce a third, very broad component 
to the Ha line. While expected to be exceedingly difficult to 
detect, its quantification will directly yield the fraction of non- 
thermal protons by energy. 

8.6. Detecting Polarization in Ha Lines 

Laming 115 711 suggests a search for linear polarization in the 
narrow Ha line, which arises when the bulk velocity of the post- 
shock flow is comparable to the thermal velocities of the elec- 
trons or protons. The polarization arising from the protons is 
always expected to be present, while the contribution from the 
electrons decreases as they heat up. 

Using the / eq parametrization (see equation fT8lO as intro- 
duced by Ghavamian et al. |32y (see their equations [1] and 
[2]), the electron and proton temperatures in a pure hydrogen 
gas are: 

3m p v 2 s l p \ 

3mrf(J_\ 
p * 16k 

The respective thermal velocities (equation H) are then: 
3v s \m p t fi \ 

3v s 1 
Vth,/? ~ — ; 

4 JTTp 

Since the bulk velocity of the post- shock fluid is Av ~ 3v s /4 
(with respect to the pre- shock gas), the thermal velocities of the 
electrons and protons are both comparable to Av when J3 <?c 1 . 
In the case of J3 = 1 , the protons still have v±, P ~ Av, but the 
electrons now have v±, e ~ 30 Av. 

The pre-shock hydrogen atoms flowing across the collision- 
less shock front "see" an anisotropic distribution of electrons 
or protons when the respective thermal velocity is ~ Av. If 
v±,e » Av, then the electronic contribution to the polarization 
is washed out by the random motions of the electrons. In Case 
A conditions, Laming estimates that polarizations of up to 10% 
may be present for BDSs with v s « 2000 km s" 1 and J3 = or 
1; the polarization in Case B conditions, for realistic J3 values, 
is expected to be under 1%. The broad Ha line is expected to 
be much less polarized, because the broad neutrals are mov- 
ing at nearly the same bulk velocity as the post-shock elec- 
trons/protons, and thus "see" nearly isotropic electron/proton 
distributions. 

If the precursor responsible for broadening the narrow lines 
(see 35.1b does not induce polarization, then the two narrow line 
contributions can be distinguished. More accurate calculations 
of the polarization present in the narrow Ha lines are contingent 
upon obtaining better cross sections for proton-hydrogen atom 
excitation (see 36.3.3b . No successful search for polarization in 
the narrow Ha lines from BDSs has ever been reported. 



8. 7. The Case for Integral Field Unit Spectroscopy 
Until now 
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spectra of BDSs have been obtained using slit 
spectroscopy or "two-dimensional" spectra — each vertical po- 
sition along the slit has its own associated spectrum. Integral 
field unit (IFU) spectroscopy yields a "data cube": one spec- 
tral and two spatial dimensions. In essence, one obtains images 
of an object at hundreds of different wavelengths (e.g., SNR 
1987A by Kjcer et al. [53]). Its strongest advantage over slit 
spectroscopy is the increased amount of spatial information. 
For example, IFU spectroscopy will be useful for studying the 
spatial variation of the widths and intensities of the narrow and 
broad Ha lines in a resolved shock filament. 
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